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Abstract 

We propose a coordinate-space regularization of the three-body prob- 
lem with zero-range potentials. We include the effective range and the 
shape parameter in the boundary condition of the zero-range potential. 
The proposed extended zero-range model is tested against atomic helium 
trimers and is shown to provide an adequate quantitative description of 
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■ these systems. 
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1 Introduction 



, The zero-range potential [Q has been extensively used over many years as a 

practically convenient form of the effective interaction. The concept employs 
separation of scales in a physical problem and allows qualitative and often quan- 
C3 ' titative description of low-energy properties of a physical system in a simple and 

^ , transparent way (see, e.g., |^). 

However, application of the zero-range potential to a three-body system 
presents a problem - a collapse of the system known as the Thomas effect 
[^. The three-body system with zero-range potentials has no ground state but 
, infinitely many bound states with vanishing spatial extension and exceedingly 

' large binding energy. 

Several attempts have been made to alleviate this problem by adding some 
sort of cutoff to the potential in momentum space ||, ||, |^ or by switching to a 
finite-range potential in certain areas of configuration space . 

We introduce an alternative coordinate space approach where the collapse is 
removed by a suitable modification of the boundary condition of the zero-range 
potential. The boundary condition is extended to include the higher order 
parameters of the effective range expansion. The three-body system acquires 
then a well defined ground state while all the simplicity and transparency of the 
zero-range model is retained. 

We apply the developed formalism to rather involved three-body systems, 
helium trimers, and show that the extended zero-range model provides an ac- 
curate description of these systems. 



2 The zero-range model and regularization 

Zero-range potentials. The quantum mechanical two-body problem with a 
zero-range potential can be formulated |l| as the free Schrodinger equation for 
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the s-state wave- function tp with the relative coordinate r and wave number k, 



(1) 



with the solution 



r?A = sin (fcr + 5{k)) , (2) 
and a boundary condition at r = expressed in terms of the scattering length 



a as 



1 d^rip) 



r%l: dr 



= k cot S{k) = — 
a 



For negative scattering length a bound state solution exists 

rip cx exp(— Kr) , 



(3) 
(4) 



where k > can be found from the boundary condition k = l/|a|- 

The zero-range model for a three-body system can be formulated as a free 
three-body wave-function ^! which satisfies the three boundary conditions 
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a|r,-rfc|* 



|r, -rfcl* a|r, -rfcl 



(5) 



|rj-rfe|=0 



where is the coordinate of i-th particle, ai is the scattering length in the 
two-body system of particles j and k with k} being a positive permutation 
of {1,2,3}. 

The derivatives in the boundary condition arc most suitably formulated 
in terms of the hyper-spheric coordinates {p,ai} (defined in the appendix): 



d 



5|r, 



rfcl 



|rj-rfc|=0 



Mi d 
p dui 



The boundary condition (||) can then be rewritten as 



1 d{a,-^) 



dai 



ai=0 



(6) 



(7) 



Hyper-spheric expansion. We shall employ the hyper-spheric adiabatic ex- 
pansion H of the three-body wave-function 



(8) 



in terms of the complete basis f2) of the solutions of the hyper-angular 

eigenvalue equation 



,2 3 



(9) 



where Vi is the potential between particles j and fc, m is the mass scale used in 
the definition of the hyper-spheric coordinates and A is the angular part|^ of the 
kinetic energy operator (see appendix). 

^Hcrc f2 is any of the tree possible sets of angles 
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For potentials without strong repulsive cores already the lowest term in the 
expansion - the so called hyper-spheric adiabatic approximation - gives a very 
good approximation to the precise solution pc| ] . Again it is the lowest term that 
causes the Thomas collapse of a three-body system with zero-range potentials. 
Therefore in the following for the sake of simplicity we shall consider only this 
problematic lowest term of the hyper-spheric expansion. Inclusion of higher 
terms is straightforward. 

The wave-function then simplifies to 

^{p.^)^^f{pmp,^), (10) 

where the hyper-radial wave- function f{p) satisfies the equation 

(-^ + ^^-Q(.)-^)/(p) = o, (11) 

where X{p) is the lowest eigenvalue in eq. £^ is the total energy and 

dm{p,n)^^p,n) . (12) 

Faddeev equations. For short-range and zero-range potentials the Faddeev 
decomposition of the angular wave- function $(/0, fl) provides a convenient frame- 
work for an analysis of the three-body system , 

i—l 

where the three components ^i{p, ai) satisfy a system of Faddeev equations 



(A - A(p)) + ^-^VMp, ^) = 0- (14) 



sin(2ai) 

Since the zero-range potentials act only on the s-waves only the latter are in- 
cluded in each of the three components fi. 

All three components of the wave-function <I> in ( p^ must be "rotated" 
into the same Jacobi system. This is done by substituting the variables and 
subsequently projecting onto the s-waves. The transformation of (pk into the 
j'th Jacobi system is given as 

1 rf-lf-'f'ik-ail 

where 



0,.=arctan + (16) 

y\j nijink J 

The expansion of ipj^k {(Xj ) for small angles aj <^ 1 reads 
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The zero-range potentials vanish identically except at the origin and we are 
therefore left with the free Faddeev equations 



52 



(18) 



which are obtained from ( p^ ) with Ix — ly = and Vt—O, and where ly^ — \ + 4:. 
The solutions are 



ipi{p,ai) ^ Ai{p) sin J^(p) (a^ - ^) 



(19) 



with the boundary condition (f>i{p^ ^) = 0. 

The factors Ai are to be determined from the boundary condition (|^) which 
can now be reformulated in terms of the angular function (f> as 



da. 



^1 



(20) 



ai=0 



The needed wave-function and its partial derivative d{ai^)/dai are 
easily obtained from (^3|) and ( p^ ) 



■■^ sm(2(pij 



,5(a,$) 



5ai 



dipi{ai 



ai=0 



+ 0(a ■) , 
".=0 ^ sin(2(/)y) ■ 



Substituting the free solutions ( |19D leads to 

AiJ/cos (^i^^^ + 



,9(a.*) 



A, 



Qi=0 



2sin[^(0,j--f)] 
sin(2(?!)y ) 



(21) 
(22) 

(23) 
(24) 



The boundary condition (20) then becomes a system of linear equations for 
the three factors A 



AiV cos 



(4)+i: 



A 



2 Sin [z. (</),, -f)] _ p 1 
sin(2(?!),y) ^/piai 



A,sin(i.|) . (25) 



A non-trivial solution exists only when the determinant of the corresponding 
matrix M(i/, p) is zero 

detM{v,p) = , (26) 

where the matrix elements are 

TTN p 1 



M, 



Ma = z^cos ^i^^^ -I- sin {1^^^ 
2 sin [i/ (^(pij — -, 



(27) 



sin(20ij) 



The solution ii{p) of the equation (^6|) provides the adiabatic potential {i^^ip) — 
1/4)1 p^ for the hyper-radial equation (11) from which one obtains the hyper- 
radial wave-function of the three-body system. 
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Asymptotic behaviour of the eigenvalues. For a system of three identical 
bosons, where ipij — tt /"i, equation ( |2^ ) simplifies to 



-z^cos(i^D 



73 



sin(z/f) 



01a 



(28) 



For large distances, p ^ a, there is a solution that asymptotically approaches 
i^(oo) = 2. Expanding (ESh in terms of 1/p around i/=2 gives the leading terms 



12 ^/Jla A 



16 3y/JIa 



(29) 



TT p p^ IT p'^ 

which is the lowest solution when no bound two-body subsystems are present. 
In this case the effective potential is of l/p^ type. 

However, when there is a two-body bound state another kind of solution 
exists for large p which asymptotically behaves as ^ ip. The leading terms 
then are 



. p 1 8 P 1 TT 
« — — r + ^ — 1= exp — ■; — r — 



A 



P 

PLO? 



p 1 16 ^ p 1 TT^ 
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The effective potential is then of Yukawa type 
A + 15/4 _ 2mB 1 



16V3 6 / Ps 
exp(-- , 

TT p 



P' 

where B = h^/(2pma^) is the two-body binding energy and b : 
The corresponding angular wave function (19) asymptotically is 



ip 



cx exp 



pa 



(30) 
(31) 

(32) 

3^|a| /tt . 
(33) 



This wave-function is non- vanishing only when a ^ ,JJi\a\ / p ^ 1. In this 
region the Jacobi coordinates x and y (defined in the appendix) are approxi- 
mately, up to the linear terms in a, given hy x ~ pa and y ~ p. This solution 
corresponds to a bound two-body state with the momentum = i/ \a\ and the 
binding energy B. The three-body wave-function then factorizes as 



(f> cx — exp(- 



■)/(2/) 



(34) 



and describes a dimer in the bound state ^ exp(— ) and a third particle with 
a relative coordinate y and the wave- function f{y). The corresponding radial 
equation asymptotically describes a two-body system with a Yukawa potential 



dp^ 



2m , 

-riE + B) 



TT p 



fip) = 



(35) 



The term — l/(4p^) in this equation cancelled the leading order term of Q{p)- 
Indeed the normalized angular Faddeev component is (asymptotically) 



(p{p, a) 



'J2p_ 



■ exp(,-p 



(36) 
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and therefore 

Q{p)^J 'P{p,a)-^ip{p,a)da = . (37) 

We have thus a correct asymptotic wave-function corresponding to a dimer and 
a third particle in a relative s-wave. 

The term Q is generally small and only is important to ensure the correct 
asymptotic behaviour. In the following practical application we shall always for 
simplicity use only the leading term — l/(4p^) instead of the full Q similar to 
the Langer correction term in Q . 

The Thomas effect and regularization. For p a the equation ( p8| ) for 
V reduces to 

TT 8 7r 

— i/cos(j/— )H r=sin(u—) = Q, (38) 

^ 2^ V3 6^ ^ ' 

which has well known imaginary roots i^o = where g = 1.006, which cause 
the Thomas and also the Efimov jlj] effects. 

These imaginary roots lead to an effective potential in the hyper-radial equa- 
tion which in the small distance region, p ^ a, is equal to {v^ — ^ 
— 1.262//9^ and the radial equation becomes 

92 1/^-1/4 2mE\ 

W^~P^ T^J 

The (negative) energy E = — ?i^K^/(2m) is negligible compared to the effective 
potential when the distance is sufhciently small, p <^ and the corresponding 
radial equation, 

92 i^g-1/4' 



f{p) = . (39) 



Qp2' p2 jfiP)-^^ (40) 

has in this region solutions of the form f{p) ^ p", where n = ^ ± i/q. For 
imaginary vq = ±ig the exponent n also acquires an imaginary part ±ig leading 
to 

f{p)oi^/pexp{±ig\np). (41) 

This wave-function has infinitely many nodes at small distances or, correspond- 
ingly, infinitely many low lying states at smaller distances. This is called the 
Thomas effect. 

A suitable modification of the boundary condition (0) is necessary in order to 
eliminate the problematic imaginary root i/q at p = which causes the Thomas 
effect. Intuitively one could generalize the zero-range potential by introducing 
the higher order terms of the effective range theory, 

, , =- + lRk^ + PR^k\ (42) 

rip ar a 2 

where R is the effective range and P is the shape parameter of the two-body 
system. This would lead to the following modification of the matrix elements 
in the eigenvalue equation ( p6| ) 

Mii = 1/ cos ^i^— ^ + sin f i^— ^ ^ 



1 diri;) 



2 



a,; 2 



(43) 
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Figure 1: The angular eigenvalue A as function of p for the ^He-trimer for 
different potential models: exponential , realistic LM2M2 and zero- 
range with P=0.13. All models have the same scattering length a=-189.05 au 
and effective range i?=13.843 au. 



and equation (j^) for three identical bosons is then replaced by an extended 
equation 



-vcos{v^) + ^sin(i/f) ^ p 



sin(z/f) 



a 2 \ p 



(44) 



This extended equation at p = has a real root v{Q) = and the Thomas 
collapse is therefore removed. 

Although the second order term with the effective range is, in principle, 
enough for elimination of the imaginary root, the fourth order term is necessary 
to ensure the correct analytic properties of the roots of the equation. 

Unlike the scattering length and effective range the parameter P has to be 
interpreted as a regularization parameter which somehow accounts for all the 
higher order terms in the expansion, rather than a true shape parameter of 
the two-body scattering. The scattering length and effective range are important 
for the correct asymptotic behaviour of the eigenvalue A at large distances while 
the P parameter accounts for the pocket region and is in this model supposed 
to absorb all the remaining short distance properties of the system. 



3 Application to helium trimers 

The helium trimer ^Hea is a challenging three-body system as there is a weakly 
bound dimer state, ^He2, where the scattering length a=-189.054 au is much 
larger than the effective range i?=13.843 au (the atomic unit of length is equal to 
the Bohr radius a_B=0. 529177 A). Numeric computations with realistic LM2M2 
potential show that there is a ground state and an extremely weakly bound 
excited state interpreted as an Efimov state (see, e.g., Q and references 
therein). 
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Figure 2: The energies of the ground and the excited state of ^He-trimer as 
function of the P parameter of the extended zero-range model. The arrows 
indicate the results of the realistic LM2M2 model ||l0| . 

For our calculations we use as in Q the mass scale m=1822.887 au (the 
atomic unit of mass is equal to the electron mass TOe=0. 510999 MeV/c^). The 
mass of the '*He atom is m(^He)=4.0026037Ti. The angular eigenvalue X{p) — 
(p) — 4 is obtained directly by numeric solution of the transcendental equation 

®- . 

On Fig. 1 we compare the angular eigenvalues obtained from the zero-range 
model with P=0.13 and from two finite-range models: the realistic LM2M2 
potential, and an exponential potential - all models having the same scattering 
length and effective range. 

At large distances, p 3> i?, the angular eigenvalues from all models approach 
each other since it is only the scattering length and effective range that deter- 
mine the asymptotic behaviour of A. 

At short distances the behaviour is different. The realistic LM2M2 model 
with a strong repulsive core produces a strongly repulsive eigenvalue. The eigen- 
value from the exponential potential converges to A(0) = as it does for all 
potentials which diverge slower than at the origin |q]. The eigenvalue from 
the zero-range model converges to A(0) = —4 according to (Q). This is precisely 
sufficient to eliminate the Thomas effect. 

Although the zero-range model gives a stronger attraction at small distances 
the correct large distance behaviour and a good overall agreement make it a solid 
alternative to finite-range potentials. One would expect that a weakly bound 
three-body state, not sensitive to the short-range details of the potential, should 
be reasonably well described by the zero-range model while one would expect 
some over-binding for strongly bound states. 

Bound states. The extended zero-range model correctly predicts the number 
of "'He-trimcr bound states - two states - for a large variation of the P param- 
eter, see Fig. 2. The energy of the weakly bound excited state, predicted rather 
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accurately, is largely independent of the P-parameter since it mostly resides in 
the outer region which is determined exclusively by the scattering length and 
effective range. The stronger bound ground state is more sensitive to the inner 
part of the effective potential, and therefore to the P-parametcr. On average the 
zero-range model gives a description similar to finite-range models, see Table 1. 

Table 1: The bound state energies of the helium trimers ^Hca and '*He2'^He 
for finite-range potential models [ p^ , and for the zero-range model with 
P=0.13. The '*He-^He scattering length is a=-189.054 au and effective range 
i?=13.843 au. For "^He-^He system a=33.261 au, i?=f8.564 au. The mass of 
■^He is m('^He)=3. 016026. For the gaussian, exponential and zero-range models 
the shown energies are calculated within the adiabatic (one-channel) approxi- 
mation (^0|). For these simple potentials the realtive accuracy of the adiabatic 
approximation is better than 1%. The LM2M2 energies are obtained with the 
full expansion (H). 



Potential 


Eoimes) (mK) 


^i^Hea) (mK) 


E^oC'Hea^He) (mK) 


LM2M2 


-125.2 


-2.269 


-13.66 


Gaussian 


-150.2 


-2.462 


-18.41 


Exponential 


-173.9 


-2.714 


-24.27 


Zero-range 


-143.7 


-2.21 


-34.0 



Radial functions. Fig. 3 compares the radial functions from different models. 
Due to the softness of the zero-range model at short distances the ground state 
wave-function is shifted to the left in comparison to the exponential and LM2M2 
potentials. However, it is still rather similar to those obtained from the finite- 
range models. The wave-functions for the excited state are roughly identical for 
all models since the spatially extended weakly bound states are less sensitive to 
the individual features of the underlying potential model. 

Non-identical particles. Another bound atomic trimer ^He2'^He is obtained 
by substituting ^He for one of the "'He atoms. The scattering length in the 
subsystem •^He-^He is a=33.261 au and the effective range i?=18.564 au. In 
order to obtain I'ip) for the system of non- identical particles we have to solve 
the general equation (^^ . 

Although the P parameter for the •^He-'^H subsystem should generally speak- 
ing be also different, we choose the same value P=0.13 for the sake of simplicity. 
With these parameters the zero-range model correctly predicts that there is only 
one bound state in this system. The binding is somewhat larger than for finite- 
range potentials but within the same range of accuracy, see Table 1. Again, we 
did not make any attempt to fit the binding energy by varying P. 

4 Conclusion 

The zero-range potential is a very useful form of effective interaction whose 
applications to three-body systems are, however, severely hampered by the 
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exponential : 




P (au) 



Figure 3: The arbitrarily normalized radial wave- functions / as function of p 
for the ground and the excited states of the ''He trimer for the potential models 
from Fig. 1. 



Thomas collapse. We propose a coordinate space regularization of the zero- 
range potential which leads to a removal of the Thomas collapse. The new 
model on one hand retains all the simplicity of the zero-range potential and on 
the other hand provides a fully regularized solution for the three-body system. 
Compared to finite-range potentials the computational load is greatly reduced 
and amounts to solving a transcendental equation for the effective potential and 
subsequent ordinary differential equation for the radial wave-function. We ap- 
ply the proposed model to atomic helium trimers and show that it works well 
and produces results comparable to finite-range models. 



A Hyper-spheric coordinates 

If TUi and Yi refer to the i-th particle then the hyper-radius p and the hyper- 



angles ai are defined in terms of the Jacobi coordinates and as [15 

mjYj + rukYk 



/i,(rj - rfe) , y, = 



1 mjUik 1 mi{mj+mk) 

M^ = , Pjk = ^ ^ (45) 

m ruj + nik m rrii + rrij + 

/9sin(aj) = , pcos(ai) = yi , 

where {i,j,k} is a cyclic permutation of {1,2,3} and m is an arbitrary mass. 
The set of angles fii consists of the hyper- angle ai and the four angles Xi/|xi| 
and yi/|yi|. The kinetic energy operator T is then given as 



sin(2Qi) daf sin^(ai) cos2(a,j) 
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where l^^^ and 1^. are the angular momentum operators related to and y,;. 
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